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In this work, we study the A a; Fe2-; / Se2 (A=K, Rb) superconductors using angle-resolved photoe- 
mission spectroscopy. In the low temperature state, we observe an orbital-dependent renormalization 
for the bands near the Fermi level in which the d xy bands are heavily renormliazed compared to the 
d X z/dy Z bands. Upon increasing temperature to above 150K, the system evolves into a state in which 
the d xy bands have diminished spectral weight while the d xz /d yz bands remain metallic. Combined 
with theoretical calculations, our observations can be consistently understood as a temperature in- 
duced crossover from a metallic state at low temperature to an orbital-selective Mott phase (OSMP) 
at high temperatures. Furthermore, the fact that the superconducting state of Aa;Fe2- H Se2 is near 
the boundary of such an OSMP constraints the system to have sufficiently strong on-site Coulomb 
interactions and Hund's coupling, and hence highlight the non-trivial role of electron correlation in 
this family of iron superconductors. 



PACS numbers: 74.25.Jb, 74.70.-b, 79.60.-i 

Electron correlation remains a central focus in the 
study of high temperature superconductors. The 
strongly correlated cuprate superconductors are under- 
stood as doped Mott insulators (MI) |1 while the iron- 
based superconductors (FeSCs) have been found to be 
at most moderately correlated [2H3] • Even though the 
low energy electronic structures of different FeSCs fam- 
ilies share the common Fe 3d bands, there are system- 
atic variations in their physical properties, such as or- 
dered magnetic moment and effective mass [5j. In gen- 
eral, electron correlation is the weakest in iron phosphides 
with relatively low mass renormalization [5] , and moder- 
ate in the more extensively studied iron arsenides [5J [3] . 
The Fe(Te,Se) chalcogenide family, in comparison, seems 
to harbor stronger correlation as inferred from larger 
ordered moment, yet metallic resistivity is still ob- 
served [BJ. The newest iron chalcogenide superconduc- 
tors, A x Fe 2 ~ v Se 2 (A=alkali metal) [7HTT] (AFS) hints at 
even stronger correlation with a large observed moment 
of 3.3/l(b |12| and insulating transport behavior in the 
phase diagram. 

Another important factor in understanding the FeSCs 
lies in their multi-orbital nature. In such a system, 
orbital-dependent behavior as well as competition be- 
tween inter- and intra-orbital interactions could play a 
critical role in determining their physical properties. One 
example is the orbital anisotropy that onsets with the 



in-plane symmetry breaking structural transition as ob- 
served in the underdoped arsenides |13j . Another exam- 
ple is the different pairing symmetry that could arise con- 
sidering the relative strength of inter- and intra-orbital 
interactions extensively studied theoretically [14]. Theo- 
retical models have considered correlation effects in the 
bad metal regime in terms of an incipient Mott pic- 
ture [5J [T5], and the proximity to the Mott transition 
may be orbital-dependent even for orbitally-independent 
Coulomb interactions p~6l419j . What arises from the 
model is an orbital selective Mott phase (OSMP), in 
which some orbitals are Mott-localized while others re- 
main itinerant. First introduced in the context of the 
Ca2_ a ;Sr 2 ;Ru04 system, an OSMP may result from both 
the orbital-dependent kinetic energy and the combined 
effects of the Hund's coupling and crystal level split- 
tings [201 HI]- An OSMP links naturally with mod- 
els of coexisting itinerant and localized electrons that 
have been proposed to compensate for the shortcom- 
ings of both strong coupling and weak coupling ap- 
proaches [221 [23]. However, to date, there has been no 
experimental evidence for OSMP in any FeSC. 

In this paper, we present angle-resolved photoemission 
spectroscopy (ARPES) data from two superconducting 
AFSs, K x Fe 2 - y Se 2 (KFS) and Rb^Fe^Sea (RFS), with 
Tc of 32K and 31K, respectively, as well as insulating 
and intermediate dopings (see SI). We observe the su- 
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FIG. 1: Measured electronic structure of K ;E Fe2-;,Se2. (a) Fermi surface mapping by integrating 20meV window about Ef- 
Green lines outline the 2-Fe Brillouin zone. (b),(d) Spectral images and (c),(e) second derivatives taken along the F-X direction 
using light polarizations as marked, (f)-(g), equivalent of that of (b)-(c) for Ba(Coo.o7Feo.93)2As2. (h) Schematic of the dominant 
orbital characters of the two electron pockets of the same size near X point, with one of the pockets (dotted) imploded for 
clarity, (i) LDA calculations [27] for KFS with the dominant orbital characters labeled, (j)-(l), Spectral images taken across 
the X-point under different polarizations and cut directions. Guides to eye for the observable bands are overlaid, with colors 
indicating the dominant orbital characters-blue: d xy ; red: d xz ; green: d yz . All data taken at 30K, with 47.5eV photons except 
(d), (e), and (k), which were taken with 26eV photons. 



perconducting AFSs undergoing a temperature-induced 
crossover from a metallic state in which all three ti g 
orbitals-d-rj,, d xz and d yz are present near the Fermi level 
(Ep) to a state in which the d xy bands has diminished 
spectral weight while the d xz /d yz bands remain metal- 
lic. In addition, the intermediate doping shows stronger 
correlation than the superconducting doping, as seen in 
the further renormalization of the d xz /d yz bands and the 
much more suppressed d xy intensity, while the insulating 
doping has no spectral weight near i?p in any orbitals. 
From comparison with our theoretical calculations, the 
ensemble of our observations are most consistent with 
the understanding that the presence of strong Coulomb 
interactions and Hund's coupling places the supercon- 
ducting AFSs near an OSMP at low temperatures, and 
crosses over into the OSMP via raised temperature, while 
the intermediate and insulating compounds are on the 
boundary of the OSMP and in the MI phase, respec- 
tively, suggesting the importance of electron correlation 
in this family of FeSCs. 

The low temperature electronic structure of 
K :r Fe2-j,Se2is shown in Fig. [I] The Fermi surface 
(FS) of KFS consists of large electron pockets at the 
Brillouin zone (BZ) corner-X point, and a small electron 
pocket at the BZ center-r point (Fig. [ija)), consistent 
with previous ARPES reports [24H26] , For the crys- 
tallographic 2-Fe unit cell, LDA calculations predict 
two electron pockets at the X point [57] (Fig. |TJi) ) . 
While FS map appears to show only one electron 
pocket at X, measurements under different polarizations 



(Fig. [Tfj)-(l)) reveal the expected two electron bands 
with nearly degenerate Fermi crossings (hp) but dif- 
ferent band bottom positions-a shallower one around 
-0.05eV and a deeper one that extends to the top of 
the d xz /d yz hole-like bands at -0.12eV. The Luttinger 
volume of the two electron pockets gives ~0.16 electrons 
per Fe. Considering that the C4 symmetry of the 
crystal dictates degeneracy of the d xz /d yz electron band 
bottom and hole band top at X, the shallower electron 
band that is not degenerate with the hole-like band at 
higher binding energy is most likely of d xy character, 
and the deeper one d xz along T-X and d yz along the 
perpendicular direction. This observed orbital character 
assignment seems to contradict the LDA prediction 
of the shallower electron band as d xz /d yz and the 
deeper one d xy in FeSC [28] (Fig. fiji)), as observed 
in Co-doped BaFe 2 As 2 (Fig. [TJf)-(g)). However, this 
assignment can be understood if we consider the KFS 
band structure as a whole. Three filled hole bands are 
observed near the T point (Figs, [ljd)-(e)), where the 
two lower ones can be identified as d xz /d yz and the 
higher one d xy . Interestingly, the d xy band is far more 
renormalized (~a factor of 10) compared to LDA than 
the d xz /d yz bands (~a factor of 3), indicating stronger 
correlation for the d xy orbital. This is in contrast to 
the Co-doped BaFe2As2 band structure, in which all 
orbitals are renormalized by roughly the same factor 
(~2) compared to LDA. Hence, our assignment of 
d xy character to the shallower electron band is consistent 
with strong orbital-dependent renormalization, which 
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FIG. 2: Temperature dependence of electron bands near X. (a) Spectral images taken with 26eV photons along cutl. (b) 
Spectral images and (c) corresponding energy distribution curves (EDCs) taken with 47.5eV photons along cut2. Polarizations 
are as marked. Spectral images for selected temperatures are shown for both geometries. To rule out a trivial thermal 
broadening for the diminishing spectral weight, we artificially introduce a 210K thermal broadening to the 30K data and plot 
them next to the 210K data for comparison. 



brings the deeper d xy electron band at X predicted 
by LDA to be shallower than the d xz /d yz band. This 
orbital-dependent renormalization behavior also emerges 
from our theoretical calculations as will be discussed 
later. 

As the electronic structure at Ep is dominated by the 
large electron pockets at X, we focus on these features. 
Fig. [2] shows a temperature-dependent study of these 
bands taken in two polarization geometries. In Fig.[2ja), 
the matrix element for d xy is stronger than d xz . At low 
temperatures, d xy band is clearly resolved with weaker 
intensity for the d xz band. With increasing temperature, 
the spectral weight of the d xy band noticeably diminishes, 
eventually revealing the remaining d xz band at high tem- 
peratures. In Fig. [2jb), the deeper d yz band has more 
intensity than the corresponding d xz band in Fig. [2^a) 
while the presence of the d xy band is still very noticeable 
from the increased intensity where they significantly over- 
lap, as well as from the energy distribution curves (EDCs) 
shown in Fig. [2jc) . With increasing temperature, again, 
the spectral weight of the d xy band diminishes, leaving 
only the d yz band at high temperatures, which clearly 
has a deeper band bottom than the shallow d xy band 
(Fig. j2jc)). In addition, we have artificially introduced a 
210K thermal broadening to the 30K spectra as shown in 
the last column of Fig. [2| The clear contrast to the 210K 
data rules out a trivial thermal broadening as an origin 
for the observed diminishing of d xy spectral weight. 

To capture this behavior quantitatively, we analyze the 
temperature dependence of the EDCs at the X-point, 
stacked in Fig. |3jc). At all temperatures, there is a large 



hump background corresponding to the large hole-like 
dispersion at high binding energy. At low temperatures, 
there is an additional peak around -0.05eV correspond- 
ing to the d xy band bottom. We fit these EDCs with a 
Gaussian background for the large hump feature and a 
Lorentzian peak for the d xy band. The integrated spec- 
tral weight for the fitted d xy peak is plotted in Fig. [3^d) , 
which decreases toward zero with increasing temperature, 
seen as a non-trivial drop between 100K and 200K. As 
an independent check against trivial thermal effect, we 
choose small regions in the raw spectral image (marked in 
Fig. |3ja)-(b)) dominated by d xy (blue), d yz (green), and 
mixed d xy and d xz /d yz (magenta/cyan) characters and 
plot their integrated intensities as a function of tempera- 
ture (Fig.[3|e)). The spectral weight from ci^-dominated 
region rapidly decreases, consistent with the fitted result 
in Fig.[3|d), while that of d yz -dominated region does not 
drop in a similar manner. The regions of mixed orbitals 
show a slower diminishing spectral weight compared to 
that for d xy , reflecting the combined contributions from 
both d xz /d yz and d xy orbitals. Although this method 
has the uncertainty of small leakage of other orbitals into 
the chosen regions, which is the cause of the finite resid- 
ual value for the d xy curve, the contrasting behavior of 
the d xy versus d xz /d yz orbitals is clearly demonstrated. 
A temperature cycle test was performed to exclude the 
possibility of sample aging (see SI) . Measurements on the 
sister compound RFS reveal similar behavior (SI), sug- 
gesting the generality in this family of superconductors. 
This observation of a selected orbital that loses coher- 
ent spectral weight while the others remain metallic is 
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FIG. 3: Quantitative analysis of the temperature dependence. (a),(b) 90K spectral image from Fig. p|a),(b). (c) EDCs at 
the X point from (a) for selected temperatures. Curves are offset for clarity. Each curve is fitted to a Gaussian (dotted line) 
for the background feature at high binding energy and a Lorentzian (shaded region) for the feature around -0.05eV. (d) The 
integrated Lorentzian spectral weight are plotted as a function of temperature, (e) Temperature-dependence of the averaged 
intensity in the regions marked by colored boxes in (a)-(b), with background (dotted box of the same energy window) for each 
box subtracted. The resulting temperature-dependent curve is then normalized by the initial value. The blue (green) region has 
dominant spectral weight from the d xy (d yz ) band, whereas the magenta (cyan) region is of mixed d xy and d xz (dyz) characters. 



reminiscent of a crossover into an OSMP in which se- 
lected orbitals become Mott localized while others remain 
metallic. 

To further understand this phenomenon, we perform 
theoretical calculations based on a five-orbital Hubbard 
model to study the metal-to-insulator transition (MIT) 
in the paramagnetic phase using a slave-spin mean-field 
method [251 OS] • At commensurate electron filling n=6 
per Fe (corresponding to Fe2 + of the parent FeSC), we 
find the ground state of the system to be a metal, an 
OSMP or a MI depending on the intra-orbital repulsion 
U and the Hund's coupling J. Furthermore, the MIT can 
be triggered by increasing temperature (Fig. |4|a)) due 
to the larger entropy of the insulating phase. At a fixed 
interaction strength (within a certain range, Fig. Qa)), 
the system goes from a metal to an OSMP and then to 
a MI with increasing temperature. The MI phase is sup- 
pressed by electron doping. By contrast, the OSMP can 
survive, as shown in Fig. [T](b) for n=6.15, which roughly 
corresponds to the filling of the superconducting state 
from ARPES measurements. From the evolution of the 
orbitally resolved quasiparticle spectral weight, Z a , as 
a function of the temperature (Fig. [4jc) ) , we show that 
the OSMP corresponds to the d xy orbital being Mott 
localized [Z — 0) and the rest of the 3d orbitals remain- 
ing delocalized (Z > 0). This result originates primarily 
from a combined effect of the orbital dependence of the 
projected density of states and the interplay between the 
Hund's coupling and crystal level splitting (see SI and 
Ref.EE])- 

To compare with the ARPES measurements, we have 
calculated the quasiparticle spectral function A(k,E) 
in the 2-Fe BZ. At low temperature and in the non- 
interacting limit U=0 (Fig.Qd)), the electronic structure 



of the model agrees well with that from LDA, with the 
d xy band deeper than the d xz /d yz band at X. This order 
switches with sufficiently strong interaction. At U=3.75 
eV (Fig. |4^e)), the d^-dominated bands are pushed 
above their d xz /d yz counterparts by the strongly orbital- 
dependent mass renormalization, as reflected in the or- 
bital dependent quasiparticle spectral weights (Fig.[4|c)). 
The mass renormalization is the largest for the d xy orbital 
(~10), and smaller for the d xz /d yz orbitals (^3), which 
is compatible with the low-temperature ARPES results 
(Fig. [lje)). The temperature- induced crossover to the 
OSMP is clearly seen from the suppression of the spectral 
weights in the d xy orbital that accompanies the reduction 
of the weights in the other orbitals (Fig. |4^c),(e)-(g)), in 
agreement with the ARPES results (Fig.T3J). 

The temperature-induced nature of the crossover con- 
strains these AFS superconductors to be very close to the 
boundary of the OSMP in the zero temperature ground 
state, which is also the superconducting state. The best 
agreement between theory and experiments is achieved at 
U~3.75 eV. While the absolute value of this interaction 
is sensitive to the parameterization of the crystal levels 
and Hund's coupling, it is instructive to make a qualita- 
tive comparison with the case of the iron pnictides. The 
enhanced correlation effects for AFS tracks the reduction 
of the width of the (U=0) d xy band, which is about 0.7 
of its counterpart in 1111 iron arsenides. 

One known concern for the AFS materials is the exis- 
tence of mesoscopic phase separation into superconduct- 
ing and insulating regions [32] . which would both con- 
tribute spectral intensity in ARPES data. From mea- 
surement on an insulating RFS sample (Fig. |S4[d)), we 
see negligible spectral weight and no well-defined disper- 
sions within O.leV from Ep, as expected for an insula- 
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FIG. 4: Theoretical calculation based on a five-orbital Hubbard model. (a),(b) Slave-spin mean-field phase diagrams of the 
five-orbital Hubbard model for KFS at J/U=0.15 and two electron densities. OSMP and MI refer to orbital-selective Mott phase 
and Mott insulator, respectively, (c) The evolution of the orbital-resolved quasiparticle spectral weight Z a with temperature at 
n=6.15, U=3.75eV, and J/U=0.15. (d)-(g) show the quasiparticle spectral functions along the F-X direction at n=6.15. The 
color curves highlight the dominant orbital of the band, with d xz /d yz in red and d xy in blue. 



tor. Hence, the insulating regions in the superconduct- 
ing compounds do not contribute spectral weight to the 
near- Ep energy range, in which the temperature-induced 
crossover is observed. We have also measured a KFS 
sample whose resistivity is intermediate between super- 



conducting and insulating (Fig. S4 b)) , and was previ- 



ously proposed to be semiconducting containing both 
metallic and insulating phases 33J . Interestingly, its 
dxzldy Z bands, which must come from the metallic phase, 
appear further renormalized by a factor of 1.3 compared 
with those of the superconducting compounds. In addi- 
tion, we resolve small but finite spectral weight for a very 



shallow d xy electron band near the X-point (Fig. S5'c)). 



As expected, the peak position is even closer to Ep, con- 
sistent with additional renormalization for the shallow 
band near X, which is harder to discern with tempera- 
ture dependent study. These observations are consistent 
with the OSMP picture in that that the metallic phase in 
this KFS compound is likely even closer to the boundary 
of OSMP at low temperatures from the mass-diverging 
behavior of the d xy bands. For the same interaction 
strength, calculation from our model also identifies the 
low temperature ground state of the superconducting, in- 
termediate, and insulating phases to be located close to 
an OSMP, just at the boundary of an OSMP, and in a 
MI phase, respectively (see SI). 

While we cannot completely rule out alternative ex- 
planations for the observations presented above, the 
consistency of the totality of the observations-including 



strongly orbital dependent band renormalization for 
d xy versus d xz /d yz at both T and X points in the low 
temperature metallic state, the non-trivial temperature- 
dependent spectral weight change for only the d xy band, 
systematic doping dependence of the related effects in the 
intermediate and insulating compounds-and the theoret- 
ical calculations makes this understanding a most likely 
scenario, suggesting that the superconductivity in this 
AFS family exists in close proximity to Mott behavior. 
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SUPPLEMENTARY INFORMATION 



Materials and methods 



Single crystals of KFS were grown by the following 
method. Precursor FeSe was synthesized using Se powder 
(Alfa 99.999%) and Fe powder (99.998%) in a 1:1 molar 
ratio. The reagents were weighed and placed in a 2mL 
alumina crucible. The quartz tube was sealed after be- 
ing flushed with argon and evacuated. The sealed quartz 
tube was placed in a furnace and heated from room tem- 



perature to 1050° C in 12 hours. The furnace remained 
at 1050°C for 12 hours, then was shut off and cooled to 
room temperature. Single crystals of KFS were obtained 
by a self- flux method with mixtures of K (Alfa, 99.95%), 
and FeSe in molar ratios of 1:3 respectively. Potassium 
has a low boiling point so a small amount of additional 
K needs to be incorporated into the growth during syn- 
thesis. The reagents were weighed and placed in a 2mL 
alumina crucible, which was then sealed in a 2mL quartz 
tube after being flushed with argon and evacuated. The 
2mL quartz tube was then placed into a larger 5ml quartz 
tube and sealed after again being flushed with argon and 
evacuated. The double sealed quartz tube technique is 
employed because potassium attacks quartz. The sealed 
quartz tubes were placed in a furnace and heated from 
room temperature to 1040° C over the course of 4 hours. 
After dwelling for 2 hours, the furnace was cooled to 
850° C. The quartz tube was then removed from the fur- 
nace, rotated 180°C, and spun in a centrifuge for a few 
seconds to separate the self-flux from the single crystals. 
Crystals with dimensions up to approximately 3mm x 
3mm x 50mm could readily be extracted from the cru- 
cible. The crystals have a platelike morphology, with 
the c-axis perpendicular to the plane. Two types of 
KFS were studied, one superconducting with TC on- 
set of 32K and chemical composition as Ko.76Fei.7 2 Se 2 , 
and the other non-superconducting with composition as 
Ko.76Fei.7sSe2. Single crystals of RFS were grown as de- 
scribed elsewhere pQ. Two types of RFS were studied: 
one superconducting with Tc onset of 31K and chemi- 
cal composition Rb .93Fei.7oSe2; and the other insulating 
with composition Rbo.goFei.7sSe2. 

As the crystals are sensitive to air, all sample prepa- 
rations were preformed inside a nitrogen-flushed glove 
box. ARPES measurements were carried out at beam- 
line 10.0.1 of the Advanced Light Source and beamlinc 
5-4 of Stanford Synchrotron Radiation Lightsource using 
SCIENTA R4000 electron analyzers. The total energy 
resolution used was 25 meV or better and the angular 
resolution was 0.3° C. Single crystals were cleaved in situ 
at low temperatures and measured in an ultra high vac- 
uum chamber with a base pressure better than 3xl0 -11 
Torr. 



Discussion of orbital assignment 

The common band structure for FeSC in the Brillouin 
zone (BZ) corresponding to the two-Fe unit cell consists 
of three hole bands near the T point and two electron 
bands at the X point [2] . Under C4 rotational symmetry, 
the d xz and d yz bands are degenerate at the in-plane 
high symmetry points of the BZ-T and X points. At 
the r point, this means that the d xz and d yz hole bands 
must be degenerate. Hence we have assigned the two 
lower degenerate hole bands in Fig. fife to d xz /d yz and the 
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higher one to d xy . At the X point, C4 symmetry dictates 
the bottom of the d xz /d yz electron band and the top of 
the d xz /dy Z hole band be degenerate. Since the shallow 
electron band is separated from the hole band by more 
than 70mcV, it is unlikely to be the d xz /d yz electron 
band, but is instead the d xy band. 

There is in principle the possibility of an alternative 
understanding. As described in Ref . [5J [5] , the alternat- 
ing arsenic atoms about the iron plane may induce parity- 
switching of certain orbitals when folding from the 1-Fe 
BZ to the 2-Fe BZ, which would swap the orbital charac- 
ters of the d xz electron band and the d xy electron band at 
X along T-X direction, making the originally d xz electron 
band to be more observable under an odd polarization 
(Fig. fib). Under this understanding, the shallow band is 
still the d xz band rather than the d xy band. However, for 
this understanding to hold, the aforementioned degener- 
acy of the d xz /d yz bands must be lifted to account for 
the 70meV gap between the shallow electron band and 
the hole band at X. 

There are two mechanisms that could lift this degen- 
eracy. One is spin-orbit coupling. However, in this com- 
pound, the magnitude of the spin-orbit coupling is not 
strong enough to cause such a large gap. Moreover, the 
effect of the spin-orbit coupling is stronger at the T than 
at the X point [3] , and should cause a larger gap between 
the d xz /d yz hole bands at T, which is not observed. The 
second mechanism is in-plane symmetry breaking, which 
causes anisotropic shift of the d xz /d yz bands and hence 
lift the degeneracy of these two orbitals, as has been ob- 
served in the orthorhombic state of BaFe2As2 [I]. Up to 
date, there has not yet been any report of in-plane sym- 
metry breaking in the AFS superconductors. In addition, 
one should note that the orthorhombicity causes the ex- 
istence of natural twinning in the crystals. Hence in an 
unstressed crystal, the ARPES signal is a combination of 
dispersions from both types of domains, and would see 
bands from T-X and T-Y directions simultaneously, as 
was observed in BaFe2As2 [I]. For the microscopically 
phase separated AFS superconductors, such twinning ef- 
fect should be even more noticeable than BaFe2As2, but 
is not observed in the dispersions of AFS superconduc- 
tors. 

Another factor that should be considered is the com- 
parison between T and X. The flat d xy hole band ob- 
served at r is renormalized by a factor of ~10 compared 
to bare LDA while the d xz /d yz hole bands are renor- 
malized by a factor of ^3. When such kind of orbital- 
dependent renormalization is considered for the bands at 
X, it would bring the d xy electron band to be shallower 
than the d xz /d yz electron bands, which is more consis- 
tent with the assignment of the shallower electron band 
to d xy . Further support for this assignment comes from 
the general consideration that a stronger mass renormal- 
ization would imply a stronger temperature dependence 
of the corresponding quasiparticle spectral weight, as has 



been shown in the main text. 

To summarize, while the possibility in principle exists 
for other unknown mechanisms and alternative explana- 
tions, given the totality of the observations and consid- 
erations given above, the understanding of the shallow 
electron band as d xy seems to be the most consistent 
overall. 



Temperature cycle test 

To test that the disappearance of the d xy orbital spec- 
tral weight with raised temperature is not due to sample 
aging, we have performed a temperature cycle test in 
which we cleaved the sample at 10K and took measure- 
ments as temperature is raised up to 300K and cooled 
back down again. Spectra images for selected tempera- 
tures as well as the EDC at X point in this temperature 
cycle are shown in Fig. |S1| where we see that the spec- 
tral weight for the d xy orbital dominated shallow elec- 
tron band diminishes at higher temperatures, and recov- 
ers when cooled back down. This shows that the observed 
temperature-induced crossover is not due to any surface 
effects and is a robust bulk phenomenon. 



Results on Rba;Fe2_ySe2 

Wc have done similar measurements on superconduct- 
ing RFS as discussed in the main text for superconduct- 
ing KFS. Fig. |S2| summarizes the main results. RFS has 
similar electronic structure as KFS, including Fermi sur- 
face, band dispersions and their orbital characters, specif- 
ically there is a shallow electron band at X that is d xy and 
a deep electron band that is d xz /d yz . Fig. S2d-f shows 
the temperature dependence of the electron bands, equiv- 
alent to Fig. [5Ja) for the KFS. Here again we see that 
the shallow d xy band is present at low temperatures, and 
its spectral weight diminishes with raised temperature. 
Similar quantitatively analysis is also done on the RFS. 



Fig. S2 'g) shows the EDC at X point taken in the geom- 
etry of Fig.fS2];d)-(f), taken from 10K to 200K. Similar 
to the case of KFS, the ^-0.05eV hump shaded by blue is 
the d xy band bottom, whose spectral weight diminishes 
with raised temperature. The EDCs are again fitted with 
a Gaussian background and a Lorentzian peak, whose 
integrated intensity as a function of temperature is plot- 
ted in Fig. S2h. The method of integrating intensity in 
boxed regions is also done on RFS. Here two regions were 
chosen (Fig. S2i): blue box for a region dominated by 
d xy band, and magenta box for a region of mixed d xy and 
d xz . The background for each region is taken as a box of 
same energy window away from dispersions, marked by 
dotted boxes in Fig. S2 d). The integrated intensity of 



these two regions as a function of temperature is plotted 
in Fig. S2 [i), showing the diminishing trend of d^region 




Ak x (u/a) -o.15 -0.10 -0.05 o.oo o.os 

E-E F (eV) 

FIG. SI: Temperature cycle test on K a; Fe2- !/ Se2 (a) Spectral 
images across X point measured at selected temperatures in 
the temperature cycle from 10K to 300K to 40K. Polarization 
geometry and photon energy is the same as that in Fig.^a). 
(b) EDC at X point taken in the temperature cycle, where the 
peak around -0.05eV is the d^band bottom, whose spectral 
weight diminishes up to 300K and recovers when cooled back 
down. 

and the slower trend of the mixed region. These behav- 
iors in RFS are very similar to those in KFS in the main 
text. 

Orbital-selective Mott phase of the five-orbital 
Hubbard model for Ka,Fe2 -ySe2 

The five-orbital Hubbard model is given by H = 
+ Hint , where H^ n and iJj nt respectively denote the 
kinetic and the on-site interaction parts of the Hamilto- 
nian. H ln t contains the intra- and inter-orbital Coulomb 
repulsion, as well as the Hund's rule coupling and the pair 
hoppings.pl)] The corresponding coupling strengths are 
respectively U, U', and J, which satisfy U' = U — 2 J. [7 
For simplicity, we consider only the density-density inter- 
actions and neglect the spin-flip and pair-hopping terms. 



The results including these terms are qualitatively the 
same. The kinetic part TJkin is a tight-binding Hamil- 
tonian, and is conveniently specified in the momentum 
space. In FeSCs, each unit cell contains two Fe ions. 
Hence, ideally the tight-binding Hamiltonian must be de- 
fined in the BZ corresponding to this two-Fe unit cell. 
However, the lattice symmetry of the FeSC system al- 
lows us to work in an unfolded BZ corresponding to one- 
Fe unit cell. Following Ref. [H] and notice that the ions 
are invariant under the transformation P Z T X and P z T yi 
where T x / y \ is the translation along x(y) direction by one 
Fe-lattice spacing, and P z refers to the reflection about 
the Fe plane, we may define a pseudocrystal momentum 
k in the extended one-Fe BZ. This pseudocrystal mo- 
mentum and the conventional momentum k are related 
by k = k + Q in d xz and d yz orbitals (where Q = (ir, 7r)), 
but k = k in other orbitals. In the extended BZ, the 
tight-binding Hamiltonian reads 

tfkin = K^(k) + (A a - M )5a/j]di a(T cW (!) 

ka/3cr 

Here £ Q ^(k) is the hopping matrix in the momentum 
space associated with orbitals a and /3, A a is the on-site 
energy reflecting the crystal field splitting, ^ is the chem- 
ical potential, and 5 a p is the Kronecker's delta function. 
The expression of £ tt( g (k) is given in Ref. [5] , and it has 
the same form as appeared in the appendix of Ref. [2J. 
We adopt the tight-binding parameters of Ref. [9] , where 
they are obtained by fitting the LDA band structure of 
KFS. To better fit the LDA results, we have further tuned 
several hopping parameters by hand from their values in 
Ref. . Using the same notation as in Ref. [5] , the tight- 
binding parameters used in this paper for KFS are listed 
in Table [I] The chemical potential corresponding to the 
electron filling n=6.15 is /j,— -0.365 eV. 

The five-orbital Hubbard model is studied using the re- 
cently developed U(l) slave-spin mean-field method [10] . 
In this method, a slave quantum S=l/2 spin is introduced 
to carry the charge degree of freedom, and the spin of the 
electron is carried by a fermionic spinon. This approach 
determines the quasiparticle spectral weight Z a in each 
orbital. An orbital a is delocalized when Z a > 0, but 
becomes Mott localized if Z a = 0. 

To compare with the ARPES data, the coherent part 
of the orbital resolved spectral function A Q (k, E) in the 
folded BZ corresponding to the two-Fe unit cell is fur- 
ther calculated by convoluting the slave-spin and spinon 
Green's functions and writing them in the Lehmann rep- 
resentation via the following formula 

^ Q (k,E)=^Z Q |Un 2 ^(E-^ k ). (2) 

Here, e^k and U^ v are respectively the v's eigenenergy 
and eigenvector of the hopping matrix ^ Q) g(k). This al- 
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FIG. S2: (a) Fermi surface measured at 10K, with integration window of 20meV about Fermi level, (b) Spectral image and 
(c) second derivative measured along F-X. Blue (green) curves trace the observed d xy (d yz ) dispersions observed under this 
polarization, (a)-(c) were measured with 47.5eV photons, (d)-(f) Spectral images measured across X in the same geometry 
as that in Fig. [2ja) for selected temperatures, (g) EDC at X measured in the same geometry as that for (d)-(f) for selected 
temperatures. Each curve is fitted to a Gaussian (dotted lines) for the background feature around -0.15eV and a Lorentzian 
(shaded region) for the feature around -0.05eV. (h) The integrated Lorentzian spectral weight is plotted as a function of 
temperature, (i) Temperature-dependence of the intensity in the regions marked by colored boxes in (d). The blue region has 
dominant spectral weight from the d xy band, whereas the magenta region is of mixed d xy and d xz . 
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TABLE I: Tight-binding parameters of the five-orbital model for K^Fe2- H Se2. Here we use the same notation as in Ref [2J. The 
orbital index «=1,2,3,4,5 correspond to d xz , d yz , d x 2_ y 2, d xy , and d 3z 2_ r 2 orbitals, respectively. The units of the parameters 
are eV. 



lows us to determine the orbital character of each band 
near the Fermi level. In Fig. |S3[ we show the spectral 
functions near the Fermi level in the five-orbital model 
at n=6.15 and T=10 K for three different U values. Com- 
pared to the U=0 band structure, the bands are strongly 



renormalized by the interactions. Moreover, with increas- 
ing U, the spectral weight of the bands with a d xy orbital 
character is reduced (Fig. S3 b)) and eventually goes 
away (Fig. |S3[ c)). The bands with the d xz /d yz char- 
acter, on the other hand, still have nonzero spectral 
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FIG. S3: Calculated quasiparticle spectral weight functions for various U. (a) Orbitally resolved quasiparticle spectral weight 
functions along the 1 to X direction of the Brillouin zone at J/U=0.15, n=6.15, T=10 K, and three different U values (in three 
columns). The total spectral weights contributed from all five orbitals are shown in the first row; the spectral weights from the 
d xy and d xz +d yz orbitals are respectively shown in the second and the third row. 



weights. This behavior clearly indicates an interaction 
driven transition to an OSMP. 



Comparison of insulating, intermediate and 
superconducting compounds 



Several factors favor stabilizing the OSMP. [T2] Firstly, 
from the orbitally-projected density of states of the non- 
interacting bands, the width of the d xy orbital is nar- 
rower than that of the other Fe 3d orbitals, especially 
the d xz /dy Z orbitals. This factor recalls the mechanism 
for the OSMP initially proposed for the Ca2_ a; Sr a ;Ru04 
system [T3]. Secondly, when the Hund's coupling is suf- 
ficiently strong compared to the (nonzero) splitting be- 
tween the d xy and d xz /d yz orbitals, the high-spin config- 
uration is favored, and, due to the crystal level splitting, 
the d xy orbital is non-degenerate and located at a higher 
energy than other orbitals. As a result, for a range of 
densities, the d xy orbital is kept at half-filling while the 
d xz /d yz orbitals are more than (though still close to) half 
filled. The non-degenerate d xy orbital has a lower repul- 
sion threshold for the Mott transition than the doubly de- 
generate d xz /d yz orbitals within a wide range of Hund's 
coupling. This picture has some connection with the one 
studied in a different regime for a three-orbital model 
away from half-filling [T3] . As a combined effect of these 
factors, the d xy orbital is more strongly localized than 
the d xz /d yz orbitals. 



Here, we compare the measured electronic structure 
of three kinds of phases in KFS and RFS. The four 
compounds shown in Fig. |S4| have compositions deter- 
mined by energy dispersive X-ray spectroscopy to be (a) 
Ko.76Fei. 72 Se2, (b) if .76Fei.7 8 Se 2 , (c) Rb . 93 Fei.7oSe2, 
and (d) Ko.9oFei.7gSe 2 . From resistivity data shown, 

(a) and (c) are superconducting, (d) is insulating, and 

(b) shows a behavior somewhat between superconducting 
and insulating, which we shall call intermediate. Its re- 
sistivity is most like the compound suggested to be semi- 
conducting in a previous report [T5] , which suggests that 
both insulating and metallic regions exist in these sam- 
ples. All compounds were measured in the same experi- 
mental geometry. The two superconducting samples ex- 
hibit similar dispersions as discussed in previous section. 
The insulating sample shows negligible spectral weight 
towards Ep and no well defined dispersions near Ep, as 
expected of an insulator. The interesting case is the inter- 
mediate sample, which has resolvable dispersions, but dif- 
ferent from that of the superconducting samples. Firstly, 
we see that the most well-resolved band is the d yz band, 
which is traced in green in the second derivative plot in 
the bottom panel. Its bandwidth is approximately renor- 
malized by a factor of 1.3 compared to its equivalent in 
the superconducting samples. Secondly, from the EDC 
taken at the X point from both the superconducting KFS 



and the intermediate KFS (Fig. S5 b)-(c)), we see that 
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K x Fe 2 . / Se 2 - superconducting K„Fe 2 . y Se 2 - intermediate q Rb„Fe 3 . y Se 2 - superconducting ^ Rb y Fe 2 . y Se 2 - insulating 




k x (rr/a) MTT/a) k x (ir/a) k„ frr/a) 



FIG. S4: Comparison of superconducting, intermediate and insulating compounds (a) Resistivity (top), ARPES spectral images 
on the F-X high symmetry cut (middle) and its second derivative plot (bottom) for superconducting Krj.76Fei.72Se2, taken at 
10K. Same measurements are shown for (b) intermediate -Ko.76Fei.7sSe2 at 10K, (c) superconducting Rbo.93Fei.7oSe2 at 10K, 
and (d) insulating Ko.9oFei.7sSe2 at 50K. All ARPES data were taken with 26eV photons in the same polarization geometry 
as that for Fig.Qd). 
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FIG. S5: Comparison of superconducting and intermediate 
KFS (a) Calculated orbital-resolved quasiparticle spectral 
weight for U=3.75eV and T=30K as a function of electron 
filling, n, for d xz /d yz (red squares) and d xy (blue triangles), 
(b) Measured EDC at X (black solid curve) from supercon- 
ducting KFS from Fig. |S4[ a), fitted to a Gaussian plus a 
Lorentzian, plotted in red and blue, respectively, (c) Same as 
(b) but for the intermediate KFS sample from Fig. S4b, with 
the EDC from the superconducting sample renormalized by 
1.3 in energy for comparison (gray dotted curve). Blue arrows 



the peak around -0.05eV indicating the d xy electron band 
bottom in the superconducting sample (Fig. [S5jb)) be- 
comes a smaller but discernible shoulder closer to Ep in 
the intermediate sample (Fig.|S5[c)). Even after a renor- 



malization of 1.3 to account for the increased renormal- 
ization of the d yz band (gray dotted line in Fig. S5 'c)) as 
discussed above, we see that the d xy band in the inter- 
mediate compound still needs a further renormalization 
going from the superconducting to the intermediate sam- 
ple (Fig.|S5|c)). Also, the relative spectral weight of the 
d xy orbital to that of d yz is much more reduced in the 
intermediate sample. Both the further renormalization 
of the bandwidth and the reduction of spectral weight 
of the d xy band indicate that the intermediate sample 
is even closer to the OSMP than the superconducting 
compounds at low temperatures. 

Assuming the same interaction strengths, the various 
phases in the KFS and RFS compounds can be under- 
stood in terms of an interplay between the electron dop- 
ing and vacancy order. In the vacancy disordered case, 
our calculation (Fig. |S5( a)) identifies a doping-induced 
transition to an OSMP near n=6.02 per Fe. Hence the 
undoped system with n=6 is already in an OSMP. The 
vacancy order further drives it through a Mott transition 
in all orbitals to a Mott insulator [16l [17] . This accounts 
for the absence of spectral weights near Ep in the in- 
sulating RFS. For the superconducting samples and the 
intermediate KFS compound, the phases giving ARPES 
signals are likely to be vacancy disordered. We therefore 
interpret them as corresponding to two vacancy disor- 
dered phases at two different electron densities (n=6.15- 
6.25 for the superconducting samples, and n=6.05-6.10 
for the metallic phase in the intermediate KFS). In addi- 
tion, we emphasize that the high temperature state of the 
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superconducting region is intrinsically different from the 
vacancy-ordered insulating phase. Rather, the supercon- 
ducting, intermediate, and insulating phases likely have 
increasing correlation as they may be located close to an 
OSMP, just at the boundary of an OSMP, and in a Mott 
insulating phase, respectively. 
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